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Abstract. A local agglomeration of cooperators can support the survival or spreading of cooperation, 
even when cooperation is predicted to die out according to the replicator equation, which is often used 
in evolutionary game theory to study the spreading and disappearance of strategies. In this paper, it 
is shown that success-driven motion can trigger such local agglomeration and may, therefore, be used 
to supplement other mechanisms supporting cooperation, like reputation or punishment. Success-driven 
motion is formulated here as a function of the game-theoretical payoffs. It can change the outcome and 
dynamics of spatial games dramatically, in particular as it causes attractive or repulsive interaction forces. 
These forces act when the spatial distributions of strategies are inhomogeneous. However, even when 
starting with homogeneous initial conditions, small perturbations can trigger large inhomogeneities by 
a pattern-formation instability, when certain conditions are fulfilled. Here, these instability conditions 
are studied for the prisoner's dilemma and the snowdrift game. Furthermore, it is demonstrated that 
asymmetrical diffusion can drive social, economic, and biological systems into the unstable regime, if these 
would be stable without diffusion. 

PACS. 02.50.Le Decision theory and game theory - 87.23.Ge Dynamics of social systems - 82.40.Ck 
Pattern formation in reactions with diffusion, flow and heat transfer - 87.23.Cc Population dynamics and 
ecological pattern formation 
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Game theory is a well-established theory of individual 
strategic interactions with applications in sociology, eco- 
nomics, and biology |1|2|3|4|5[6] . and with many publica- 
tions even in physics (see Ref. [7] for an overview). It dis- 
tinguishes different behaviors, so-called strategies i, and 
expresses the interactions of individuals in terms of pay- 
offs Pij. The value Pij quantifies the result of an interac- 
tion between strategies i and j for the individual pursuing 
strategy i. The more favorable the outcome of the inter- 
action, the higher is the payoff Pij. 

There are many different games, depending on the 
structure of the payoffs, the social interaction network, 
the number of interaction partners, the frequency of in- 
teraction, and so on |4l5j . Theoretical predictions for the 
selection of strategies mostly assume a rational choice ap- 
proach, i.e. a payoff maximization by the individuals, al- 
though experimental studies |8l9ll0j support conditional 
cooperativity [TT] and show that moral sentiments [T2] can 
support cooperation. Some models also take into account 
learning (see, e.g. |13j and references therein), where it 
is common to assume that more successful behaviors are 
imitated (copied). Based on a suitable specification of the 
imitation rules, it can be shown [14115116] that the result- 
ing dynamics can be described by game-dynamical equa- 



tions |17U8j , which agree with replicator equations for the 
fitness-dependent reproduction of individuals in biology 
|19I2()I21] , 

Another field where the quantification of human be- 
havior in terms of mathematical models has been ex- 
tremely successful concerns the dynamics of pedestrians 
[2"2"] , crowds [25] , and traffic [M] . The related studies have 
led to fundamental insights into observed self-organization 
phenomena such as stop-and-go waves [24] or lanes of 
uniform walking direction [25]. In the meantime, there 
are many empirical [25] and experimental results [2 6127] . 
which made it possible to come up with well calibrated 
models of human motion [28129] . 

Therefore, it would be interesting to know what hap- 
pens if game theoretical models are combined with models 
of driven motion. Would we also observe self-organization 
phenonomena in space and time? This is the main question 
addressed in this paper. Under keywords such as "assort- 
ment" and "network reciprocity" , it has been discussed 
that the clustering of cooperators can amplify coopera- 
tion, in particular in the prisoner's dilemma [30 31 32 33 . 
Therefore, the pattern formation instability is of prime 
importance to understand the emergence of cooperation 
between individuals. In Sec. [4] we will study the in- 
stability conditions for the prisoner's dilemma and the 
snowdrift game. Moreover, we will see that games with 
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success-driven motion and asymmetrical diffusion may 
show pattern formation, where a homogeneous distribu- 
tion of strategies would be stable without the presence 
of diffusion. It is quite surprising that sources of noise 
like diffusion can support the self-organization in systems, 
which can be described by game-dynamical equations with 
success-driven motion. This includes social, economic, and 
biological systems. 

We now proceed as follows: In Sec. [2] we introduce 
the game-dynamical replicator equation for the prisoner's 
dilemma (PD) and the snow-drift game (SD). In partic- 
ularly, we discuss the stationary solutions and their sta- 
bility, with the conclusion that cooperation is expected 
to disappear in the prisoner's dilemma. In Sec. [3j we ex- 
tend the game-dynamical equation by the consideration 
of spatial interactions, success-driven motion, and diffu- 
sion. Section 3.1 compares the resulting equations with 



reaction-diffusion-advection equations and discusses the 
similarities and differences with Turing instabilities and 
differential flow-induced chemical instabilities (DIFICI). 
Afterwards, Sec. [4] analyzes the pattern formation insta- 
bility for the prisoner's dilemma and the snowdrift game 
with its interesting implications, while details of the insta- 
bility analysis are provided in Appendix [B] Finally, Sec. 
[5] studies the driving forces of the dynamics in cases of 
large deviations from stationary and homogeneous strat- 
egy distributions, before Sec. [6] summarizes the paper and 
presents an outlook. 



2 The Prisoner's Dilemma without Spatial 
Interactions 

In order to grasp the major impact of success-driven mo- 
tion and diffusion on the dynamics of games (see Sec. [3]), 
it is useful to investigate first the game-dynamical equa- 
tions without spatial interactions. For this, we represent 
the proportion of individuals using a strategy i at time t by 
Pi(t). While the discussion can be extended to any num- 
ber of strategies, we will focus on the case of two strategies 
only for the sake of analytical tractability. Here, i = 1 shall 
correspond to the cooperative strategy, i = 2 to defection 
(cheating or free- riding) . According to the definition of 
probabilities, we have < Pi(t) < 1 for i e {1,2} and the 
normalization condition 



Pi{t)+p 2 (t) = 1. 



(1) 



Let Pij be the payoff, if strategy i meets strategy j. Then, 
the expected payoff for someone applying strategy i is 



£?<(*) = £P iiPj -(t), 



(2) 



as pj (t) represents the proportion of strategy j, with which 
the payoffs must be weighted. The average payoff in 
the population of individuals is 

2 22 

E(t) = ]r m) P i(t) = ^^K(i)%-(t) . (3) 

(=1 1=1 3 = 1 



In the game-dynamical equations, the temporal increase 
dpi (t) j dt of the proportion of individuals using strategy 
i is proportional to the number of individuals pursuing 
strategy i who may imitate, i.e. basically to Pi(t). The 
proportionality factor, i.e. the growth rate \(i,t), is given 
by the difference between the expected payoff Ei{t) and 
the average payoff E(t): 



dpijt) 
dt 



\(i,t) Pi (t) = [Ei(t) - E{t)] Pi (t) 



3 = 1 1=1 j=l ' 

The equations (|4 1 are known as replicator equations. They 
were originally developed in evolutionary biology to de- 
scribe the spreading of "fitter" individuals through their 
higher reproductive success |19I20I21] , However, the repli- 
cator equations were also used in game theory, where 
they are called "game-dynamical equations" |17ll8j . For 
a long time, it was not clear whether or why these equa- 
tions could be applied to the frequency pi (t) of behavioral 
strategies, but it has been shown that the equations can be 
derived from Boltzmann-like equations for imitative pair 
interactions of individuals, if "proportional imitation" or 
similar imitation rules are assumed [14115116] . 

Note that one may add a mutation term to the right- 
hand side of the game-dynamical equations Q . This term 
could, for example, be specified as 

W 2 p 2 (t) - W lPl (t) = rq[l-p x {t)\ - r(l - q)pi{t) 

= r[q-p 1 {t)] (5) 

for strategy i = 1, and by the negative expression of this 
for i = 2. Here, r is the overall mutation rate, W 2 = rq the 
mutation rate towards cooperation, and W± = r(l — q) the 
mutation rate towards defection |16j . This implementation 
reflects spontaneous, random strategy choices due to erro- 
neous or exploration behavior and modifies the stationary 
solutions. 

It can be shown that 



E 

i=l 



dpijt) 
dt 



(6) 



so that the normalization condition (JT|) is fulfilled at all 
times t, if it is fulfilled at t = 0. Moreover, the equa- 
tion dpi(t)/dt = \piit) implies Pi(t) > at all times t, if 
Pi(0) > for all strategies i. 

We may now insert the payoffs of the prisoner's 
dilemma, i.e. 



Pn = R ("reward"), 
P12 = S ("sucker's payoff"), 
P21 = T ("temptation"), and 
P 22 = P ("payoff") 



with the assumed payoff relationships 
T > R> P > S . 



(7) 



(8) 



Dirk Helbing: Pattern Formation in Games with Success-Driven Motion 



3 



Additionally, one often requires 

2R> S + T. 



(9) 



The "reward" R is the payoff for mutual cooperation 
and the "punishment" P the payoff for mutual defection, 
while T is the "temptation" of unilateral defection, and a 
cheated cooperator receives the sucker's payoff S. While 
we have P > S in the prisoner's dilemma, the snowdrift 
game (also known as chicken or hawk-dove game) is char- 
acterized by S > P, i.e. it is defined by 



T > R> S > P . 



(10) 



Both games are characterized by a temptation to defect 
(T > R), while the prisoner's dilemma has the additional 
challenge that there is a high risk to cooperate (S < P). 
This difference has a large influence on the resulting level 
Pi(t) of cooperation: Inserting the above payoffs into Eq. 
B|, one eventually obtains the game-dynamical equation 



d P i(t) 
dt 



= [1 - Pl (t)}[ - A + B Pl (t)] Pl (t) . 



(11) 



This directly follows from Eq. ( 46 1 of Appendix [A] when 
the abbreviations 



A= P - S and B=P+R-S-T 



(12) 



are used to pronounce the equation's structure. Setting 
dpi(t)/dt = 0, one obviously finds three stationary solu- 
tions p\ (t) — p\ , namely 



0. 



p\ = \ and p\ 



A 
B 



(13) 



Not all of these solutions are stable with respect to small 
deviations. In fact, a linear stability analysis (see Ap- 
pendix [A| delivers the eigenvalues 



Ai 



-A, 



A, 



A-B, 



and 



A :! 



A 1 



A' 

B, 
(14) 

The stationary solution p\ is stable with respect to per- 
turbations (i.e. small deviations from them), if Afc < 0, 
while for > 0, the deviation will grow in time. Due to 



A- B = T - R> 0, 



(15) 



the solution p\ = 1 corresponding to 100% cooperators 
is unstable with respect to perturbations, i.e. it will not 
persist. Moreover, the solution p\ will be stable in the 
prisoner's dilemma because of ^4 = P — 5 > 0. This cor- 
responds to 0% cooperators and 100% defectors, which 
agrees with the expected result for the one-shot prisoner's 
dilemma (if individuals decide according to rational choice 
theory). In the snowdrift game, however, the stationary 
solution p\ = is unstable due to A = P — S < 0, while 
the additional stationary solution p\ = A/B < 1 is sta- 
ble. Hence, in the snowdrift game with B < A < 0, we 
expect the establishment of a fraction A/B of coopera- 
tors. For the prisoner's dilemma, the solution p\ does not 



exist, as it does not fall into the range between and 
1 that is required from probabilities: If B > 0, we have 
p\ = A/B > 1, while p\ = A/B < for B < 0. 

In summary, for the prisoner's dilemma, there is no 
evolutionarily stable solution with a finite percentage of 
cooperators, if we do not consider spontaneous strategy 
mutations (and neglect the effect of spatial correlations 
through the applied factorization assumption). According 
to the above, cooperation in the PD is essentially expected 
to disappear. Strategy mutations, of course, can increase 
the stable level p\ of cooperation from zero to a finite 
value. Specifically, p\ will assume a value close to zero for 
small values of r, while it will converge to q in the limit 
r — > oo. 

In the next section, we will show that 

1. when the proportions of cooperators and defectors are 
allowed to vary in space, i.e. if the distribution of coop- 
erators and defectors is inhomogeneous, the proportion 
of cooperators may locally grow, 

2. we obtain such a variation in space by success-driven 
motion, as it can destabilize a homogeneous distribu- 
tion of strategies, which gives rise to spatial pattern 
formation in the population (agglomeration or segre- 
gation or both [49]). 

Together with the well-known fact that a clustering of co- 
operators can promote cooperation |30l31l32l33j . pattern 
formation can potentially amplify the level of cooperation, 
as was demonstrated numerically for a somewhat related 
model in Ref. [53][j]We will also show that, in contrast 
to success-driven motion, random motion ("diffusion" in 
space) stabilizes the stationary solution p\ with 0% coop- 
erators (or, in the presence of strategy mutations, with a 
small percentage of cooperators) . 



3 Taking into Account Success- Driven 
Motion and Diffusion in Space 

We now assume that individuals are distributed over dif- 
ferent locations x £ [0, L] of a one-dimensional space. A 
generalization to multi-dimensional spaces is easily possi- 
ble. For the sake of simplicity, we assume that the spatial 
variable x is scaled by the spatial extension L, so that x is 
dimensionless and varies between and 1. In the follow- 
ing, the proportion of individuals using strategy i at time 
t and at a location between x and x + dx is represented 
by pi(x,t)dx with pi(x,t) > 0. Due to the spatial degrees 



1 In contrast to this EPJB paper, the one published in Ad- 
vances in Complex Systems (ACS) studies the dynamics in 
a two-dimensional grid, assuming spatial exclusion (i.e. a cell 
can only be occupied once), neglecting effects of noise and dif- 
fusion, and choosing the payoffs P = S = 0, which restricts the 
results to a degenerate case of the prisoner's dilemma and the 
snowdrift game. Moreover, this EPJB paper focusses on the 
pattern formation instability rather than the amplification of 
the level of cooperation, formalizes social forces resulting from 
success-driven motion, and discusses pattern formation in the 
spatial prisoner's dilemma induced by asymmetrical diffusion. 
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of freedom, the proportion of defectors is not immediately 
given by the the proportion of cooperators anymore, and 
the previous normalization condition p 2 (t) = 1 — pi (t) is 
replaced by the less restrictive condition 



2 

E 



dxpi(x, t) 



(16) 



This allows the fractions of cooperators and defectors to 
uncouple locally, i.e. the proportion of cooperators does 
not have to decrease anymore by the same amount as the 
proportion of defectors increases. 

Note that, if pi(x,t) = pi(x,t)N/L represents the den- 
sity of indi vidu als pursuing strategy i at location x and 
time i, Eq. ( 16 1 can be transferred into the form 



2 

E 

i=l 



dx Pi(x, t) 



N 
1 



(17) 



where N is the total number of individuals in the system 
and p their average density. 

One may also consider to treat unoccupied space for- 
mally like a third strategy i = 0. In this case, however, 
the probabilities pi(x,t) in all locations x add up to the 
maximum concentration p max , see Eq. (23 1. This means 



and 



Po(x,t) = p max -pi(x,t) -p 2 (x,t) 

dp (x,t) ^ ggipM) _ dp 2 (x,t) 
dt ~ " at &t 



(18) 



(19) 



Therefore, po(x 7 t) can be eliminated from the system of 
equations, because unoccupied space does not interact 
with strategies 1 and 2. As a consequence, the dynamics 
in spatial games with success-driven motion is different 
from the cyclic dynamics in games considering volunteer- 
ing |35j : In order to survive invasion attempts by defectors 
in the prisoner's dilemma, one could think that coopera- 
tors would seek separated locations, where they would be 
"loners" . However, cooperators do not tend to maneuver 
themselves into non-interactive states [S3]: On the con- 
trary: The survival of cooperators rather requires to have 
a larger average number of interaction partners than de- 
fectors have. 

After this introductory discussion, let us now extend 
the game-dynamical equations according to 



dpi(x,t) 
dt 



2 2 2 2 v 

E pi E Pi 3 p j ~ EE^^j W x ' *) 

1=1 3 = 1 1=1 3 = 1 ' 



D 



dx\ M ' } dx 

with the local expected success 



Ei(x,t) = y^ i P l0 p j {x,t) 

3 = 1 



d 2 p l (x,t) 
dx 2 



(20) 



(21) 



compare Eq. (J2|. The additional sum X^pz had to be in- 
troduced for reasons of normalization, as we do not have 
Pi + P2 = 1 any longerj^] dpi( x, t)/dt represents the (par- 
tial) time derivative. T he fi rst term in large brackets on the 
right-hand side of Eq. ( pp| assumes that locally, an imita- 
tion of more successful strategies occurs. An extension of 
the model to interactions with neighboring locations would 
be easily possible. The second term, which depends on 
Ei(x,t), describes success-driven motion |36l37j . Finally, 
the last term represents diffusion, and Di > are called 
diffusion coefficients or diffusivities. These terms can be 
generalized to multi-dimensional spaces by replacing the 
spatial derivative d/dx by the nabla operator V. 

The notion of success-driven motion is justified for the 
following reason: Comparing the term describing success- 
driven motion with a Fokkcr-Planck equation [38 , one can 
conclude that it corresponds to a systematic drift with 
speed 

dEj(x,t) 

Vi{x,t) = — . (22) 

According to this, individuals move into the direction of 
the gradient of the expected payoff, i.e. the direction of 
the (greatest) increase of Ei(x, t). In order to take into ac- 
count capacity constraints (saturation effects), one could 
introduce a prefactor 



C(z,i) = l-E 



Pi(x,t)N 



i=i 



Pn 



>o, 



(23) 



where p max = p m a,xN/ L > N / L > represents the maxi- 
mum density of individuals. This would have to be done 
in the imitation-based replicator terms in the first line of 
Eq. (20 1 as well. In the following, however, we will focus 
on the case C = 1, which allows for a local accumulation 
of individuals. 

The last term in Eq. ( 20 ) is a diffusion term which re- 
flects effects of random motion in space [38] . It can be eas- 
ily seen that, for Di > 0, the diffusion term has a smooth- 
ing effect: It eventually reduces the proportion pi(x,t) in 
places x where the second spatial derivative d 2 pi/dx 2 is 
negative, in particular in places x where the distribution 
Pi(x,t) has maxima in space. In contrast, the proportion 
Pi(x,t) increases in time, where d 2 pi/dx 2 > 0, e.g. where 
the distribution has its minima. Assuming an additional 
smoothing term 

d 4 pi(x,t) 



D 



dx* 



(24) 



2 Rather than multiplying the first sum over j in Eq. ( 20 1 by 
one could also divide the second sum over j by 2j ; Pi> 
corresponding to a subtraction of the average expected suc- 
cess ~^2i Pi Ei/ ~^2i Pi from the expected success Ei. The alter- 
native specification chosen here assumes that the number of 
strategic game-theoretical interactions of an individual per unit 
time is proportional to the number of individuals it may inter- 
act with, i.e. proportional to ^2,pi(x, t). Both specifications 
are consistent with the game-dynamical equation Q, where 
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with a small constant Do > on the right-hand side of 
Eq. ( 20 ) makes the numerical solution of this model well- 



behaved (see Appendix |B|. 

When Eq. (20) is solved, one may, for example, assume 
periodic boundary conditions (i.e. a circular space). In 
this case, we have Pi(l,t) = Pi(0,t) and d k p i (l 1 t)/dx k = 
d k pi(0, t) /dx k , and by means of partial integration, it can 
be shown that 







i 1 

qI^Z J dxpi(x,t) = 0. 
i=i 



(25) 



Therefore, Eq. ( 20 I fulfils the normalization condition ( 16 1 
at all times, if it is satisfied at t = 0. Furthermore, it can 
be shown that Pi(x, t) > for all times t, if this is true at 
time t = for all strategies i and locations x. 



3.1 Comparison with Reaction- Diffusion- Ad vection 
Equations 

It is noteworthy that the extended game-dynamical model 
( 20 ) has some similarity with reaction-diffusion-advection 



(RDA) equations. These equations have been developed to 
describe the dynamics of chemical reactions with spatial 
gradients, considering the effects of differential flows and 
diffusion. Specifically, the kinetics of (binary) chemical re- 
actions is reflected by non-linear terms similar to those in 
the first line on the right-hand side of Eq. ( p0| ) . The first 
term in the second line represents ad vection terms (dif- 
ferential flows) , while the last term of Eq. ( 20 I delineates 
diffusion effects. 

Considering this apparent similarity, what dynamics 
do we expect? It is known that reaction-diffusion equa- 
tions can show a Turing instability [39 40 (also without 
an advection term). Specifically, for a chemical activator- 
inhibitor system, one can find a linearly unstable dynam- 
ics, if the diffusivities D\ and D 2 are different. As a con- 
sequence, the concentration of chemicals in space will be 
non-homogeneous. This effect has been used to explain 
pattern formation processes in morphogenesis |41l42j . Be- 
sides the Turing instability, a second pattern-forming in- 
stability can occur when chemical reactions are coupled 
with differential flows. These so-called "differential flow- 
induced chemical instabilities" (DIFICI) can occur even 
for equal or vanishing diffusivities Di = D, but they may 
also interact with the Turing instability [43 4 4145146 137] . 

The difference of the extended game-dynamical equa- 
tion (20 1 as compared to the RDA equations lies in the 



specification of the velocity Vi(x,t), which is determined 
by the gradient of the expected success Ei(x,t). There- 
fore, the advective term is self-generated by the success 
of the players. We may rewrite the corresponding term in 
Eq. pOl as follows: 



d_ 

Ox 



dx 



^ PlJ dx 



dpi dpj 
dx 



13 dx 2 



According to the last term of this equation, the advec- 
tion term related to success-driven motion implies effects 
similar to a diffusion term with negative diffusion coef- 
ficients PiPij (if Pij > 0). Additionally, there is a non- 
linear dependence on the gradients dpi/dx and dpj/dx of 
the strategy distributions in space. Both terms couple the 
dynamics of different strategies j. 

As a consequence of this, the resulting instability con- 
ditions and dynamics are different from the RDA equa- 
tions. We will show that, without success-driven mo- 
tion, diffusion cannot trigger pattern formation. Diffusion 
rather counteracts spatial inhomogeneities. Success-driven 
motion, in contrast, causes pattern formation in a large 
area of the parameter space of payoffs, partly because of 
a negative diffusion effect. 

Besides this, we will show that different diffusivities Di 
may trigger a pattern formation instability in case of pay- 
off parameters, for which success-driven motion does not 
destabilize homogeneous strategy distributions in the ab- 
sence of diffusion. This counter-intuitive effect reminds of 
the Turing instability, although the underlying mathemat- 
ical model is different, as pointed out before. In particular, 
the largest growth rate does not occur for finite wave num- 
bers, as Appendix[B] shows. The instability for asymmetric 
diffusion is rather related to the problem of "noise-induced 
transitions" in systems with multiplicative noise, which 
are characterized by space-dependent diffusion coefficients 
[4"B] . A further interesting aspect of success-driven motion 
is the circumstance that the space-dependent diffusion ef- 
fects go back to binary interactions, as the multiplicative 
dependence on pi and pj shows. 



4 Linear Instability of the Prisoner's Dilemma 
and the Snowdrift Game 

After inserting the payoffs ^ of the prisoner's dilemma 
and the snowdrift game into Eq. ( |20| , one can see the 
favorable effect on the spreading of cooperation that the 
spatial depend ence , in particular the relaxe d n ormaliza- 
tion condition (16) can have: For i = l, Eq. (20) becomes 



dp!(x, t) 
dt 



= P2 

d_ 

dx 



(R - T) Pl + 
{D 1 ~p l R) 



(S - P) P 2 

dpi 



dx 



PiS 



Pi 

dp 2 
dx 



(27) 



It is now an interesting question, whether the agglomer- 
ation of cooperators can be supported by success-driven 
motion. In fact, the second line of Eq. (27 1 can be rewrit- 
ten as 



(£>i - Pl R) 



d 2 Pi 
dx 2 



R 



dpi 
dx 



PiS 



d 2 P2 c,dpi dp 2 



dx 2 



S 



dx 



(26) 



dx 

This shows that a curvature d 2 pi / dx 2 < can support the 
increase of the proportion of strategy i as compared to the 
game-dynamical equation Q without spatial dependence. 
The situation becomes even clearer, if a linear stability 
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analysis of Eqs. (27 1 is performed (see Appendix IbI . The 



result is as follows: If the square of the wave number k, 
which relates to the curvature of the strategy distribution, 
is large enough (i.e. if the related cluster size is sufficiently 
small), the replicator terms in the first line on the right- 
hand side of equation (20 1 become negligible. Therefore, 
the conditions, under which homogeneous initial strategy 
distributions are linearly unstable, simplify to 



and 



p\R + p 2 P >Di+D 2 



(plR-D^P-D,) Kp^ST, 



(29) 



(30) 



see Eqs. ^ andj75])0 



If condition (I29I) or condition (30) is fulfilled, we ex- 



pect emergent spatio-temporal pattern formation, basi- 
cally agglomeration or segregation or both [49 . As has 
been pointed out before, the agglomeration of cooperators 
can increase the level of cooperation. This effect is not pos- 
sible, when spatial interactions are neglected. In that case, 
we stay with Eq. (Ill, and the favorable pattern-formation 



effect cannot occur. 

Let us now discuss a variety of different cases: 

1. In case of diffusive motion only (i.e. no success-driven 
motion) , Eq. ( 29 I must be replaced by 



> Di + D 2 



(31) 



and Eq. (30 I by 



D r D 2 <0, (32) 

which cannot be fulfilled. Therefore, diffusion without 
success-driven motion does not support pattern forma- 
tion or a related increase in the level of cooperation. 
In the case of success-driven motion with finite diffu- 



sion, Eqs. (29 1 and (30 I imply 



P ° 1 R+(l-p° 1 )P> D 1 +D 2 



(33) 



and 



(p?J?-A)[(l-rf)P-£> 2 ] <p° 1 (l-p° 1 )ST, (34) 

if the normalization condition p 2 = 1 — pi for a 
homogeneous initial condition is taken into account. 
These instability conditions hold for both, the pris- 
oner's dilemma and the snowdrift game. Inequality 
(33 1 basically says that, in order to find spontaneous 



pattern formation, the agglomerative tendency R of co- 
operators or the agglomerative tendency P of defectors 
(or both) must be larger than the diffusive tendency. 
This agglomeration, of course, requires the reward R 



This instability condition has been studied in the context 
of success-driven motion without imitation and for games with 
symmetrical payoff matrices (i.e. Pij = Pji), which show a 
particular behavior [49]. Over here, in contrast, we investigate 
continuous spatial games involving imitation (selection of more 
successful strategies), and focus on asymmetrical games such 
as the prisoner's dilemma and the snowdrift game (see Sec.|4|, 
which behave very differently. 



of cooperation to be positive, otherwise cooperators 
would not like to stay in the same location. The al- 
ternative instability condition ( 34 I requires that the 



product RP of the payoffs resulting when individuals 
of the same strategy meet each other is smaller than 
the product ST of payoffs resulting when individuals 
with different strategies meet each other. This basi- 
cally excludes a coexistence of the two strategies in the 
same location and is expected to cause segregation. It 
is noteworthy that condition ( 34 1 is not invariant with 

by a constant value 



respect to shifts of all payoffs P^ 
c, in contrast to the replicator equations ( 46 ) 



In case of the prisoner's dilemma without strategy 
mutations, the stable stationary solution for the case 
without spontaneous strategy changes is pi = p\ = 0. 
This simplifies the instability conditions further, yield- 
ing 



and 



P> D 1 + D 2 



D 1 (P-D 2 )<0. 



(35) 



(36) 



In order to fulfil one of these conditions, the punish- 
ment P must be positive and larger than D 2 to support 
pattern formation (here: an aggregation of defectors). 
Naturally, the survival or spreading of cooperators re- 
quires the initial existence of a finite proportion pi > 
of them, see the previous case. 

In the special case P = S = [34 , Eqs. ((33} and ([Ml 
become 

piR >D 1 + D 2 (37) 



and 



(p° 1 R-D 1 )D 2 <0, 



(38) 



which requires piR > D\. Therefore, a finite initial 
proportion pi of cooperators is needed again for pat- 
tern formation (an agglomeration of cooperators) . This 
can be easily reached by spontaneous strategy changes. 
5. Neglecting diffusion for a moment (i.e. setting D\ = 
D 2 = 0), no pattern formation should occur, if the 
condition 

RP > ST (39) 



and, at the same time, 



p1r + {i-p1)p <0 



(40) 



is fulfilled. Equation (39 1 implies the stability condi- 
tion 

S<f. (41) 

Besides T > R, we have to consider here that S < P 
in the prisoner's dilemma and S > P in the snowdrift 
game (see Fig.nJF] 



^Strictly speaking, we also need to take into account Eq. 
( 40 1, which implies P < for the stationary solution p? = p\ — 



of the prisoner's dilemma and generally P < —p{R/(l —p\). 
In case of the snowdrift game |50| . it is adequate to insert the 
stationary solution p\ = p\ = A/B, which is stable in case 
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Finally, let us assume that the stability conditions 
ST-RP < and p\R+ (1 -p?)P < for the previous 
case without diffusion (Di = D 2 = 0) are fulfilled, so 
that no patterns will emerge. Then, depending on the 
parameter values, the instability condition following 
fromEq.pl, 



-{D 1 -D 2 ){l-p\)P 

>o 

< D 2 Jp 1 R+(l-p 1 )P] 

>o <o 

-D^ +p^l - p\) {ST - RP) , 

<o 



(42) 



>o 



<o 



may still be matched, if (D\ — D 2 ) is sufficiently large. 
Therefore, aymmetrical diffusion (Di D 2 ) can trig- 
ger a pattern formation instability, where the spatio- 
temporal strategy distribution without diffusion would 
be stable. The situation is clearly different for symmet- 
rical diffusion with D\ = D 2 , which cannot support 
pattern formation. 

Although the instability due to asymmetrical diffusion 
reminds of the Turing instability, it must be distin- 
guished from it (see Sec. |3.1[ ). So, how can the insta- 
bility then be explained: The reason for it may be 
imagined as follows: In order to survive and spread, co- 
operators need to be able to agglomerate locally and 
to invade new locations. While the first requirement 
is supported by success-driven motion, the last one is 
promoted by a larger diffusivity D\ > D 2 of coopera- 
tors. 



5 Social Forces in Spatial Games with 
Success-Driven Motion 

In the previous section, we have shown how success-driven 
motion destabilizes homogeneous strategy distributions in 
space. This analysis was based on the study of linear 
(in)stability (see Appendix |b|) . But what happens, when 
the deviation from the homogeneous strategy distribution 
is large, i.e. the gradients dpi/dx are not negligible any 
longer? This can be answered by writing Eq. ( 22 ) explic- 
itly, which becomes 



Vi(x,t) 



3=1 



dpj{x,t) 
dx 



(x,t). 



(43) 



of no spatial interactions. This leads to the condition P < 
-plR/(l-pl) = AR/(A-B) = (P — S)R/(T — R), i.e. SR < 
(2R — T)P. The question is, whether this condition will reduce 
the previously determined area of stability given by S < RP/T 
with P < 0, see Eq. (41 1 and Fig. [I] This would be the case, 
if (2 - T/R)P > RP/T, which by multiplication with RT/P 
becomes (2RT-T 2 ) > R 2 or (T-Rf < 0. Since this condition 
cannot be fulfilled, it does not impose any further restrictions 
on the stability area in the snowdrift game. 



s 



Snow- 
drift 

game (SD) 




Prisoner's 

Dilemma 

(PD) 



S=P 



S=RP/T 



no pattern formation 

Fig. 1. Payoff-dependence of pattern-formation in the pris- 
oner's dilemma with S < P and the snowdrift game with 
S > P according to a linear stability analysis for spatial games 
with success-driven motion, but no diffusion, strategy muta- 
tions, or noise. One can clearly see that spontaneous pattern- 
formation prevails (green area), and that there is only a small 
area for P < (marked red), where a homogeneous initial 
condition is stable with respect to small perturbations. 



Here, the expression 



fij {x, t) Pi 



dpj(x,t) 
' dx 



(44) 



(which can be extended by saturation effects), may be in- 
terpreted as interaction force ("social force") excerted by 
individuals using strategy j on an individual using strat- 
egy zj^J It is visible that the sign of P\j determines the 
character of the force. The force is attractive for posi- 
tive payoffs Pij > and repulsive for negative payoffs 
Pij < 0. The direction of the force, however, is determined 
by spatial changes dpj(x,t)/dx in the strategy distribu- 
tion pj(x,i) (i.e. not by the strategy distribution itself). 

It is not the size of the payoffs P^ which determines 
the strength of the interaction force, but the payoff times 
the gradient of the distribution of the strategy j one in- 
teracts with (and the availability and reachability of more 
favorable neighboring locations, if the saturation prefactor 



5 Note that this identification of a speed with a force is some- 
times used for dissipative motion of the kind m a d 2 x a /dt 2 — 
—'yadxa/dt + ^2 g F a f3(t), where x a (t) is the location of an in- 
dividual a, the 'inass" m a reflects inertia, j a is a friction co- 
efficient, and F a f3(t) are interaction forces. In the limiting case 
m a — > 0, we can make the adiabatic approximation dx a /dt = 
YjpFaPit)/!* = X^/^M' where dx a /dt is a speed and 
fupit) are proportional to the interaction forces F a is(t). Hence, 
the quantities fcp (t) are sometimes called "forces" themselves. 
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C is taken into account). Due to the dependence on the 
gradient dpj(x,t)/dx, the impact of a dispersed strategy 
j on individuals using strategy i is negligible. This partic- 
ularly applies to scenarios with negative self-interactions 
(P i3 -<0). 

Note that success-driven motion may be caused by re- 
pulsion away from the current location or by attraction 
towards more favorable neighborhoods. In the prisoner's 
dilemma, for example, cooperators and defectors feel a 
strong attraction towards areas with a higher proportion 
of cooperators. However, cooperators seek each other mu- 
tually, while the attraction between defectors and cooper- 
ators is weaker. This is due to T + S < 2R, see inequality 
(J9j. As a result, even if P > 0, cooperators are moving 
away from defectors due to R > P in order to find more 
cooperative locations, while defectors are following them. 

Another interesting case is the game with the payoffs 
P\i = P-22 = —P and P12 = P21 = Q > P, where we have 
negative self-interactions among identical strategies and 
positive interactions between different strategies. Simula- 
tions for the no-imitation case show that, despite of the 
dispersive tendency of each strategy, strategies tend to ag- 
glomerate in certain locations thanks to the stronger at- 
tractive interactions between different strategies (see Fig. 
3 in Ref. [37]). 

The idea of social forces is long-standing. Montroll 
used the term to explain logistic growth laws 51J, and 
Lewin introduced the concept of social fields to the so- 
cial sciences in analogy to electrical fields in physics [52] . 
However, a formalization of a widely applicable social force 
concept was missing for a long time. In the meantime, so- 
cial forces were successfully used to describe the dynamics 
of interacting vehicles [23] or pedestrians [22], but there, 
the attractive or repulsive nature was just assumed. At- 
tempts to systematically derive social forces from an un- 
derlying decision mechanism were based on direct pair in- 
teractions in behavioral spaces (e.g. opinion spaces), with 
the observation that imitative interactions or the readi- 
ness for compromises had attrative effects (1 61153] . Here, 
for the first time, we present a formulation of social forces 
in game-theoretical terms. Considering the great variety of 
different games, depending on the respective specification 
of the payoffs Py, this is expected to find a wide range 
of applications, in particular as success-driven motion has 
been found to produce interesting and relevant pattern 
formation phenomena [37134] . 

6 Summary, Discussion, and Outlook 

In this paper, we have started from the game-dynamical 
equations (replicator equation) , which can be derived from 
imitative pair interactions between individuals |14ll5j . It 
has been shown that no cooperation is expected in the 
prisoner's dilemma, if no spontaneous strategy mutations 
are taken into account, otherwise there will be a signifi- 
cant, but usually low level of cooperation. In the snowdrift 
game, in contrast, the stationary solution corresponding 
to no cooperation is unstable, and there is a stable solution 
with a finite level of cooperation. 



These considerations have been carried out to illus- 
trate the major difference that the introduction of spa- 
tial interactions based on success-driven motion and dif- 
fusion makes. While diffusion itself tends to support ho- 
mogeneous strategy distributions rather than pattern for- 
mation, success-driven motion implies an unstable spatio- 
temporal dynamics under a wide range of conditions. As 
a consequence, small fluctations can destabilize a homo- 
geneous distribution of strategies. Under such conditions, 
the formation of emergent patterns (agglomeration, segre- 
gation, or both) is expected. The resulting dynamics may 
be understood in terms of social forces, which have been 
formulated here in game-theoretical terms. 

The destabilization of homogeneous strategy distribu- 
tions and the related occurence of spontaneous pattern 
formation has, for example, a great importance for the 
survival and spreading of cooperators in the prisoner's 
dilemma. While this has been studied numerically in the 
past [34] , future work based on the model of this paper and 
extensions of it shall analytically study conditions for the 
promotion of cooperation. For example, it will be interest- 
ing to investigate, how relevant the imitation of strategies 
in neighboring locations is, how important is the rate of 
strategy changes as compared to location changes, and 
how crucial is a territorial effect (i.e. a limitation p max of 
the local concentration of individuals, which may protect 
cooperators from invasion by defectors). 

Of course, instead of studying the continuous game- 
dynamical model with success-driven motion and calculat- 
ing its instability conditions, one can also perform agent- 
based simulations for a discretized version of the model. 
For the case without imitation of superior strategies and 
symmetrical payoffs (Pjj = Pji), it has been shown that 
the analytical instability conditions surprisingly well pre- 
dict the parameter areas of the agent-based model, where 
pattern- formation takes place 49J. Despite the difference 
in the previously studied model (see footnote 2), this is 
also expected to be true for the non-symmetrical games 
studied here, in particular as we found that the influence 
of imitation on the instability condition is negligible, if 
the wave number k characterizing inhomogeneities in the 
initial distribution is large. This simplified the stability 
analysis a lot. Moreover, it was shown that asymmetri- 
cal diffusion can drive our game-theoretical model with 
success-driven motion from the stable regime into the un- 
stable regime. While this reminds of the Turing insta- 
bility [39], it is actually different from it: Compared to 
reaction-diffusion-advection equations, the equations un- 
derlying the game-dynamical model with success-driven 
motion belong to another mathematical class, as is elabo- 
rated in Sec. 13.11 

In Sec. [5] it was pointed out that, in the prisoner's 
dilemma, cooperators evade defectors, who seek cooper- 
ators. Therefore, some effects of success-driven motion 
(leaving unfavorable neighborhoods) may be interpreted 
as punishment of the previous interaction partners, who 
are left behind with a lower overall payoff. However, 
"movement punishment" of defectors by leaving unfavor- 
able environments is different from the "costly" or "altru- 
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istic punishment" discussed in the literature [54]: In the 
strict sense, success-driven motion neither imposes costs 
on a moving individual nor on the previous interaction 
partners. If we would introduce a cost of movement, it 
would have to be paid by both, cooperators who evade de- 
fectors, and defectors who follow them. Therefore, costly 
motion would be expected to yield similar results as be- 
fore, but it would still be different from altruistic punish- 
ment. It should also be pronounced that, besides avoid- 
ing unfavorable locations, success-driven motion implies 
the seeking of favorable environments, which has nothing 
to do with punishment. Without this element, e.g. when 
individuals leave unfavorable locations based on a ran- 
dom, diffusive motion, success-driven motion is not effec- 
tive in promoting cooperation. Therefore, the mechanism 
of success-driven motion, despite some similar features, is 
clearly to be distinguished from the mechanism of punish- 
ment. 

Finally, note that migration may be considered as one 
realization of success-driven motion. Before, the statistics 
of migration behavior was modeled by the gravity law 
[55156] or entropy approaches |57l58j . while its dynamics 
was described by partial differential equations |59|60J and 
models from statistical physics [61 . The particular poten- 
tial of the approach proposed in this paper lies in the inte- 
gration of migration into a game-theoretical framework, as 
we formalize success-driven motion in terms of payoffs and 
strategy distributions in space and time. Such integrated 
approaches are needed in the social sciences to allow for 
consistent interpretations of empirical findings within a 
unified framework. 
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A Linear Stability Analysis of the 

Game- Dynamical Equation Without Spatial 

Interactions 

Inserting the payoffs ([7| of the prisoner's dilemma or the 
snowdrift game into Eq. Q, we get the game-dynamical 
equation 



dpijt) 
alt 



[R Pl + Sp 2 - R( Pl ) 2 

(S + T) Pl p 2 ~P(p 2 ) 2 ] Pl (t). (45) 



Considering Eq. (|T|), i.e. p 2 (t) = 1 — Pi{t), we find 
- = [Rp 1 + S(l-p 1 )-R( Pl ) 2 



dp 
~dt 



-(S + T) Pl (l - Pl ) - P(l - Pl ) 2 ] Pl (t) 
= {1- P1 )[(S-P) + (P + R-S- T) Pl ] Pl (t) . (46) 

This is a mean- value equation, which assumes a factoriza- 
tion of joint probabilities, i.e. it neglects correlations [16 . 
Nevertheless, the following analysis is suited to provide 



insights into the dynamics of the prisoner's dilemma and 
the snowdrift game. Introducing the useful abbreviations 



A = P-S 



and 



B = P + R-S-T, (47) 



Eq. (46 1 can be further simplified, and we get 

.- Pl (t)][-A + B Pl (t)] Pl (t). (48) 



dpi(t) 



dt 

Obviously, shifting all payoffs Pjj by a constant value c 
does not change Eq. (JTTJ) , in contrast to the case involving 
spatial interactions discussed later. 

Let p\ with k S {1,2,3} denote the stationary solu- 
tions ( 13 1 of Eq. (Ill, defined by the requirement dpi /dt = 
0. In order to analyze the stability of these solutions with 
respect to small deviations 



Spi(t)= Pl (t) 



Pl ■ 



(49) 



we perform a linear stability analysis in the following. For 
this, we insert Eq. ( 49 ) into Eq. Ill, which yields 



dSpijt) 
dt 



= {l-p k -5p x ){-A + Bp k +B5p x ){p k 

= [(l-p k )p k + (l-2p k )S Pl -(Sp 1 ) 2 ] 
x (-A + Bp k +B8px). 



Spi) 



(50) 



If we concentrate on sufficiently small deviations 5p± (t) , 
terms containing factors [Spi(t)] m with an integer expo- 
nent m > 1 can be considered much smaller than terms 
containing a factor 8pi{i). Therefore, we may linearize the 
above equations by dropping higher-order terms propor- 
tionally to [5pi(t)] m with m > 1. This gives 

= [(1 - p k )p\ + (1 - 2p k )S Pl (t)} {-A + Bp k ) 



+(l-p k )p k B6 Pl (t) 
[(l-2p k )(-A + Bp k ) 



(l-p k )p k B]6 Pl (t),(51) 



as (1 — p k )p\{— A + Bp k ) = for all stationary solutions 
p\. With the abbreviation 



A, = (1 - 2p k )(-A + Bp k ) + (1 - p\)p\B . 



we can write 



dS Pl (t) 
dt 



(52) 



(53) 



If Afc < 0, the deviation 8p\{t) will exponentially decay 
with time, i.e. the solution will converge to the stationary 
solution p k , which implies its stability. If A^, > 0, however, 
the deviation will grow in time, and the stationary solution 
p\ is unstable. For the statio nary solutions p\ = 0, p\ = 1, 
and Pi — A/B given in Eq. ( 13 1, we can easily find 



A, : . 1. A, = i B . and A, = . 1 ( 1 - - ) . 



respectively. 



(54) 
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B Linear Stability Analysis of the Model with 
Success-Driven Motion and Diffusion 

In order to understand spatio-temporal pattern formation, 
it is not enough to formulate the (social) interaction forces 
determining the motion of individuals. We also need to 
grasp, why spatial patterns can emerge from small per- 
turbations, even if the initial distribution of strategies is 
uniform (homogeneous) in space. 

If Eq. ( p0| is written explicitly for i — 1, we get 

dpi 
dt 



ipl +P2)(PllPl + Pl2P2) - P\l{Pl) 2 
-(Pl2 + P2l)piP2 ~ P 22 (p 2 f]pi(x,t) 



d 

dx 



Pi [Pi 



dpi 
dx 



Pi 



(¥2 
' dx 



cPpi 
dx 2 



(55) 



The equation for i — 2 looks identical, if only pi{x,t) 
and p2{x,t) are exchanged in all places, and the same is 
done with the indices 1 and 2. We will now assume a 
homogeneous initial condition pi(x,0) — pf (i.e. a uni- 
form distribution of strategies i in space) and study the 
spatio-temporal evolution of the deviations Spi(x,t) = 
Pi(x,t) — pf. Let us insert for p\ one of the values p\, 
which are sta tionary solutions of the partial differential 
equation (20 1, as p\ = (1 — pi) holds for homogeneous 



stra tegy distributions due to the normalization condition 
(16). Assuming small deviations Spi(x,t) and linearizing 
Eq. (55 1 by neglecting non- linear terms, we obtain 



dSpi 



(Pll - P2l)plp 2 + (Pl2 - P22){pl) 2 
(Pn-P2l)(p° 1 SP2+P°2Sp 1 ) 

-2(P 12 -P 22 )p° 2 Sp 2 
d 2 5p 



Spi 



Pi 



P°i P11- 



dx 2 



Pi 



d 2 Sp 2 
dx 2 



Di 



d 2 S Pl 
dx 2 



(56) 



Again, a mutation term W 2 Sp2(x, t) — WiSp\(x, t) reflect- 
ing spontaneous strategy changes may be added, see Eq. 
d5f. The analogous equation for Sp 2 (x,t) is obtained by 
exchanging strategies 1 and 2. 

In Eq. ( 56 1 , it can be easily seen that success-driven 
> has a similar functional form, but 



motion wit 



Pi 



the opposite sign as the diffusion term. While the latter 
causes a homogenization in space, success-driven motion 
can cause local agglomeration [37], first of all for Py > 0. 

It is known that linear partial differential equations 
like Eq. (56 1 are solved by (a superposition of) functions 
of the kina 



5pi(x,t) 



ai cos(kx) + hi sin(K.T) 



(57) 



where and bi arc initial amplitudes, A — A(k) is their 
growth rate (if positive, or a decay rate, if negative), and 
k = K n = 2nn/L with n £ {1,2, . . .} are possible "wave 
numbers" . The "wave length" 2tt/k = L/n may be imag- 
ined as the extension of a cluster of strategy i in space. 



Obviously, possible wave lengths in case of a circular space 
of diameter L — 1 are fractions L/n. The general solution 
of Eq. ([56]) is 



5pi(x, t) = } j e A(Kre) * a i>n cos{n n x) + b ijn s'm(K n x) 

71=1 

_ (58) 

i.e. a linear superposition of solutions of the form (57) with 
all possible wave numbers k„. For t = 0, t he exponential 
prefactor e A ( K ")* becomes 1, and Eq. (58 1 may then be 



viewed as the Fourier series of the spatial dependence of 
the initial condition Sp(i, x, 0). Hence, the amplitudes Oi jn 
and bi n correspond to the Euler- Fourier coefficients [r?2"] . 

Let us now de termine the possible eigenvalues A(k). 
For the ansatz (57l, we have dSpi(x,t)/dt — XSpi(x,t) and 
d 2 Spi(x,t) /dx 2 — —K 2 Spi(x,t). Therefore, the linearized 
equations can be cast into the following form of an eigen- 
value problem with eigenvalues A: 



A 



'spi(x,ty 

Sp 2 (x,t) 



M21 M 22 



5pi(x,t) 
Sp 2 (x,t) 



(59) 



=M 



Here, we have introduced the abbreviations 
M n = A u + (p°P u - D^n 2 , 
Mu - A 12 + pIPuk 2 , 

M 2 1 = A21 + p° 2 P2lK 2 , 

M 22 = A 22 + (p° 2 P22 - D 2 )K 2 

with 



An = 


(P12 


Au = 


(P11 


A21 = 


(P22 


A22 = 


(P21 



Pi 



Pi 



p% 



Pi 



(60) 
(61) 
(62) 
(63) 

(64) 
(65) 
(66) 
(67) 



The eigenvalue problem (59 1 can only be solved, if the 



determinant of the matrix (M — Al) vanishes, where 1 de- 
notes the unit matrix |62j . In other words, A are solutions 
of the so-called "characteristic polynomial" 

(M n - A)(M 22 - A) - M 12 M 21 
= A 2 - (Mn + M 22 )A + M n M 22 - M 12 M 21 = . (68) 

This polynomial is of degree 2 in A and has the following 
two solutions: 



A(«) 



Mn + M 22 



± ^V(Mn + M22) 2 - 4(MnM 22 - M 12 M 2 i) . (69) 

The fastest growing mode (i.e. the value of the wave num- 
ber k with the largest real value Re(A) of A) usually de- 
termines the length scale of the emerging patterns. Con- 
sidering ( 60 1 to ( 63 1 , we can easily see that the largest 
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value of Re (A) is reached in the limit k — > oo. This is due 
to the relationship of success-driven motion with negative 
diffusion. Hence, the finally resulting distribution would 
be a superposition of delta peaks. As this is not favor- 
able from a numerical perspective, the smoothing term 



(24) may be added, which implies the additional terms 



Dqk 4 in Eqs. (60 1 and (63 1. These terms imply, in fact 
that Re(A) reaches its maximum value for a finite value of 
k. Note, however, that discrete models involving success- 
driven motion also tend to end up with distributions ap- 
proximating a superposition of delta peaks |36|37j 



When deriving the instability conditions from Eq. ( 69 I 
in the following, we will focus on the particularly interest- 
ing case, where the mathematical expression under the 
root is non-negative (but the case of a negative value if 
4M 12 M 2 i < -(Mn - M 22 ) 2 could, of course, be treated 
as well). It can be shown that A becomes positive, if one 
of the following instability conditions is fulfilled: 

M n + M 22 > (70) 

or 

MnM 22 < M 12 M 21 . (71) 

In this case, we expect the amplitudes of the small de- 
viations 5pi(x,t) to grow over time, which gives rise to 
spatial pattern formation (such as segregation). Inserting 



the abbreviations (60 1 to (63 1, the instability conditions 
become 

[A 11 + {jP x P 11 -D 1 )k 2 ] + [A 22 + {p q 2 P 22 -D 2 )k 2 ] >0 (72) 
and 

[An + (p?Pn - D^k 2 ] [A 22 + (p° 2 P 22 - D 2 )k 2 ] 
< {A 12 + p\Pi 2 n 2 ) (A 21 + p° 2 P 21 k, 2 ) . (73) 

If k is large enough (i.e. if the related cluster size is suf- 



ficiently small) , the instability conditions ( 72 I and ( 73 I 
simplify to 

p\Pii+plP 22 >D 1 + D 2 (74) 

and 

(rfPu - £>i)(^P 22 - D 2 ) < P \P 12 p° 2 P 21 . (75) 



These are further discussed in the main text. 



